Translating Answer-Set Programs into Bit- Vector Logic* 



< 



Mai Nguyen, Tomi Janhunen, and Ilkka Niemela 

Aalto University School of Science 
Department of Information and Computer Science 

{Mai . Nguyen, Tomi . Janhunen, Ilkka . Ni erne 1 a} @ aalto . f i 



(._:> . Abstract. Answer set programming (ASP) is a paradigm for declarative problem solving where prob- 

ed ' lems are first formalized as rule sets, i.e., answer-set programs, in a uniform way and then solved by 

j^jQ, computing answer sets for programs. The satisfiability modulo theories (SMT) framework follows a 

1-^ ' similar modelling philosophy but the syntax is based on extensions of propositional logic rather than 

rules. Quite recently, a translation from answer-set programs into difference logic was provided — 
enabling the use of particular SMT solvers for the computation of answer sets. In this paper, the trans- 
' lation is revised for another SMT fragment, namely that based on fixed-width bit-vector theories. Thus, 

even further SMT solvers can be harnessed for the task of computing answer sets. The results of a 
preliminary experimental comparison are also reported. They suggest a level of performance which is 
similar to that achieved via difference logic. 

Q ; 1 Introduction 

Answer set programming (ASP) is a rule-based approach to declarative problem solving [15,22,24]. The 
^ . idea is to first formalize a given problem as a set of rules also called an answer-set program so that the 

' answer sets of the program correspond to the solution of the problem. Such problem descriptions are 

CO i typically devised in a uniform way which distinguishes general principles and constraints of the problem in 

00 ' question from any instance-specific data. To this end, term variables are deployed for the sake of compact 

. representation of rules. Solutions themselves can then be found out by grounding the rules of the answer- 

OO ' set program, and by computing answer sets for the resulting ground program using an answer set solver 

. State-of-the-art answer set solvers are already very efficient search engines [7, 11] and have a wide range 

of industrial applications. 

The satisfiability modulo theories (SMT) framework [3] follows a similar modelling philosophy but 
the syntax is based on extensions of propositional logic rather than rules with term variables. The SMT 
. framework enriches traditional satisfiability (SAT) checking [5] in terms of background theories which are 

' selected amongst a number of alternatives.' Parallel to propositional atoms, also theory atoms involving 

^ non-Boolean variables^ can be used as references to potentially infinite domains. Theory atoms are typically 

used to express various constraints such as linear constraints, difference constraints, etc., and they enable 
very concise representations of certain problem domains for which plain Boolean logic would be more 
verbose or insufficient in the first place. 

As regards the relationship of ASP and SMT, it was quite recently shown [20,25] that answer-set 
programs can be efficiently translated into a simple SMT fragment, namely difference logic (DL) [26]. 
This fragment is based on theory atoms of the form x — y < k formalizing an upper bound k on the 
difference of two integer-domain variables x and y. Although the required transformation is linear, it is not 
reasonable to expect that such theories are directly written by humans in order to express the essentials of 
ASP in SMT. The translations from [20, 25] and their implementation called LP2diff'' enable the use of 
particular SMT solvers for the computation of answer sets. Our experimental results [20] indicate that the 
performance obtained in this way is surprisingly close to that of state-of-the-art answer set solvers. The 
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results of the third ASP competition [7], however, suggest that the performance gap has grown since the 
previous competition. To address this trend, our current and future agendas include a number of points: 

- We gradually increase the number of supported SMT fragments which enables the use of further SMT 
solvers for the task of computing answer sets. 

- We continue the development of new translation techniques from ASP to SMT. 

- We submit ASP-based benchmark sets to future SMT competitions (SMT-COMPs) to foster the effi- 
ciency of SMT solvers on problems that are relevant for ASP. 

- We develop new integrated languages that combine features of ASP and SMT, and aim at implementa- 
tions via translation into pure SMT as initiated in [18]. 

This paper contributes to the first item by devising a translation from answer-set programs into theories 
of bit-vector logic. There is a great interest to develop efficient solvers for this particular SMT fragment 
due to its industrial relevance. In view of the second item, we generalize an existing translation from [20] to 
the case of bit-vector logic. Using an implementation of the new translation, viz. LP2bv, new benchmark 
classes can be created to support the third item on our agenda. Finally, the translation also creates new 
potential for language integration. In the long run, rule-based languages and, in particular, the modern 
grounders exploited in ASP can provide valuable machinery for the generation of SMT theories in analogy 
to answer-set programs: The source code of an SMT theory can be compacted using rules and term variables 
[18] and specified in a uniform way which is independent of any concrete problem instances. Analogous 
approaches [2, 14, 23] combine ASP and constraint progranaming techniques without a translation. 

The rest of this paper is organized as follows. First, the basic definitions and concepts of answer-set pro- 
grams and fixed-width bit- vector logic are briefly reviewed in Section 2. The new translation from answer- 
set programs into bit- vector theories is then devised in Section 3. The extended rule types of smodels 
compatible systems are addressed in Section 4. Such extensions can be covered either by native transla- 
tions into bit-vector logic or translations into normal programs. As part of this research, we carried out a 
number of experiments using benchmarks from the second ASP competition [11] and two state-of-the-art 
SMT solvers, viz. BOOLECTOR and z3. The results of the experiments are reported in Section 5. Finally, 
we conclude this paper in Section 6 in terms of discussions of results and future work. 

2 Preliminaries 

The goal of this section is to briefly review the source and target formalisms for the new translation devised 
in the sequel. First, in Section 2.1, we recall normal logic programs subject to answer set semantics and the 
main notions exploited in their translation. A formal account of bit- vector logic follows in Section 2.2. 

2.1 Normal Logic Programs 

As usual, we define a normal logic program P as a finite set of rules of the form 

o-s- 6i,...,6„,~ci,...,~Cm (1) 

where a, &i, . . . , 6„, and ci, . . . , Cm are propositional atoms and ~ denotes default negation. The head of 
a rule r of the form (1) is hd(r) = a whereas the part after the symbol forms the body of r, denoted 
by bd(r). The body bd(r) consists of the positive part bd'''(r) = {bi, . . . , 6„} and the negative part 
bd~(r) = {ci, . . . , Cm} so that bd(r) = bd~'~(r) U {~c | c G bd~(r)}. Intuitively, a rule r of the form (1) 
appearing in a program P is used as follows: the head hd(r) can be inferred by r if the positive body atoms 
in bd'''(r) are inferable by the other rules of P, but not the negative body atoms in bd~ (r). The positive 
part of the rule, r"*" is defined as hd(r) ■(— bd'''(r). A normal logic program is called positive if r = r+ 
holds for every rule r € P. 



Semantics To define the semantics of a normal program P, we let At(P) stand for the set of atoms that 
appear in P. An interpretation of P is any subset / C At(P) such that for an atom a € At(P), a is true in 
I, denoted I \= a, iff a G I. For any negative literal ^c, / |= iff J ^ c iff c ^ I. A rule r is satisfied in 
/, denoted / |= r, iff / |= bd(r) implies / \= hd(r). An interpretation / is a classical model of P, denoted 
7 1= P, iff, 7 ^ r holds for every r e P. A model M |= P is a minimal model of P iff there is no M' \= P 
such that M' C M. Each positive normal program P has a unique minimal model, i.e., the least model 
of P denoted by LM(P) in the sequel. The least model semantics can be extended for an arbitrary normal 
program P by reducing P into a positive program P^ = {r+ | r e P and M n bd~ (r) =0} with respect 
to M C At(P). Then answer sets, also known as stable models [16], can be defined. 

Definition 1 (Gelfond and Lifschitz [16]). An interpretation M C At(P) is an answer set of a normal 

program P iff M = LM(P*=^). 

Example 1. Consider a normal program P [20] consisting of the following six rules: 

a 6, c. a i— d. b i— a, 

6 a, ~c. c d ~c. 

The answer sets of P are Mi = {a, b, d} and M2 = {c}. To verify the latter, we note that P^^ = {a ^ 
b,c; b <^ a; c a ^ d} for which LM(P^=) = {c}. On the other hand, we have P'^^ = P^^ for 
M3 = {a, 6, c} so that M3 ^ AS(P). ■ 

The number of answer sets possessed by a normal program P can vary in general. The set of answer sets 
of a normal program P is denoted by AS(P). Next we present some concepts and results that are relevant 
in order to capture answer sets in terms of propositional logic and its extensions in the SMT framework. 

Completion Given a normal program P and an atom a G At(P), the definition of a in P is the set of rules 
Defp(a) = {r G P I hd(r) = a}. The completion of a normal program P, denoted by Comp(P), is a 
propositional theory [8] which contains 

a^\J{/\bA/\^c) (2) 

reDefp(a) 6ebd+(r) cebd-(r) 

for each atom a G At(P). Given a propositional theory T and its signature At(T), the semantics of T is 
determined by CM(T) = {M C At(T) | M ^ T}. It is possible to relate CM(Comp(P)) with the models 
of a normal program P by distinguishing supported models [1] for P. A model M |= P is a supported 
model of P iff for every atom a E M there is a rule r £ P such that hd(r) = a and M [= bd(r). In general, 
the set of supported models SuppM(P) of a normal program P coincides with CM(Comp(P)). It can be 
shown [21] that stable models are also supported models but not necessarily vice versa. This means that in 
order to capture AS(P) using Comp(P), the latter has to be extended in terms of additional constraints as 
done, e.g., in [17,20]. 

Example 2. For the program P of Example 1, the theory Comp(P) has formulas a •(-> {b t\c)\/ d,b ^ 
(a A -irf) V (a A -ic), c -f-^ -irf, and d ^ -ic. The models of Comp(P), i.e., its supported models, are 
Ml = {a, b, d}, M2 = {c}, and M3 = {a, b, c). ■ 

Dependency Graphs The positive dependency graph of a normal program P, denoted by DG"'" (P), is a pair 
(At(P), <) where b < a holds iff there is a rule r € P such that hd(r) = a and b G bd'''(r). Let <* denote 
the reflexive and transitive closure of <. A strongly connected component (SCC) of DG^(P) is a maximal 
non-empty subset 5 C At(P) such that a <* b and b <* a hold for each a,b € S. The set of defining 
rules is generalized for an SCC S by Defp(S') — Uaes Dcfp(a). This set can be naturally partitioned into 
sets Extp(5) = {r G Defp(S') | bd+(r) n 5 = 0} and Intp(S') = {r G Defp(S') | bd+(r) n 5 9^ 0} of 
external and internal rules associated with S, respectively. Thus, Def p{S) = Ext p (5) U Int p {S) holds in 
general. 



Example 3. In the case of the program P from Example 1, the SCCs of DG~''(P) are Si = {a, b}, S2 = 
{c}, and 5*3 = {d}. For Si, we have Extp(S'i) = {o d}. ■ 



2.2 Bit- Vector Logic 



Fixed-width bit-vector theories have been introduced for high-level reasoning about digital circuitry and 
computer programs in the SMT framework [27, 4J. Such theories are expressed in an extension of proposi- 
tional logic where atomic formulas speak about bit vectors in terms of a rich variety of operators. 

Syntax As usual in the context of SMT, variables are realized as constants that have a free interpretation 
over a particular domain (such as integers or reals)''. In the case of fixed-width bit-vector theories, this 
means that each constant symbol x represents a vector x[l . . .m] of bits of particular width m, denoted 
by w(a;) in the sequel. Such vectors enable a more compact representation of structures like registers 
and often allow more efficient reasoning about them. A special notation n is introduced to denote a bit 
vector that equals to n, i.e., n provides a binary representation of n. We assume that the actual width 
m > log2 (n + 1) is determined by the context where the notation n is used. For the purposes of this paper, 
the most interesting arithmetic operator for combining bit vectors is the addition of two m-bit vectors, 
denoted by the parameterized function symbol +m in an infix notation. The resulting vector is also m-bit 
which can lead to an overflow if the sum exceeds 2™ — 1. Moreover, we use Boolean operators =m and 
<m with the usual meanings for comparing the values of two m-bit vectors. Thus, assuming that x and y 
are m-bit free constants, we may write atomic formulas Uke x y and x <m y in order to compare the 
m-bit values of x and y. In addition to syntactic elements mentioned so far, we can use the primitives of 
propositional logic to build more complex well-formed formulas of bit- vector logic. The syntax defined for 
the SMT hbrary contains further primitives which are skipped in this paper. A theory T in bit- vector logic 
is a set of well-formed bit- vector formulas as illustrated by the following example. 

Example 4. Consider a system of two processes, say A and B, and a theory T = {a {x <2 y), h 
{y <2 x)} formalizing a scheduling pohcy for them. The intuitive reading of a (resp. b) is that process A 
(resp. B) is scheduled with a higher priority and, thus, should start earlier. The constants x and y denote the 
respective starting times of A and B. Thus, e.g., x <2 y means that process A starts before process B. ■ 

Semantics Given a bit-vector theory T, we write At(r) and FC(T) for the sets of propositional atoms and 
free constants, respectively, appearing in T. To determine the semantics of T, we define interpretations for 
T as pairs (/, t) where / C At(T) is a standard propositional interpretation and r is a partial function 
that maps a free constant x € FC(T) and an index 1 < i < w{x) to the set of bits {0, 1}. Given r, a 
constant x G FC(T) is mapped onto t{x) = '^'ili\T{x, i) ■ 2^(^)^*) and, in particular, r(n) = n for any 
n. To cover any well-formed terms^ ti and t2 involving +m and m-bit constants from FC(r), we define 
T{ti -\-m ^2) = T{ti) + T{t2) mod 2™ and w(ti +,„ ^2) = ni. Hence, the value T{t) can be determined 
for any weU-formed term t which enables the evaluation of more complex formulas as formalized below. 

Definition 2. Let T be a bit-vector theory, a G At(T) a propositional atom, ti and t2 well-formed terms 
over FC(T) such that w(<:i) = w{t2), and (j) and ip well-formed formulas. Given an interpretation {I,t) 
for the theory T, we define 

1. {I, r) 1= a <S=4> a e I, 

2. {I,t) \= ti =m t2 ^ T(ii) = T{t2), 

3. {I,t) h h <m t2 ^ T{ti) < T{t2), 

4. {I,t)^^<P ^ {I,r)^<l>, 

5. {I,T)^4>y^ ^ {I,T)\=cl>or{I,T)^^, 

6. {I,t) [=(j)^tp ^=> {I, t) l^cjjor (J, t) \= tP, and 

7. {I,t) \= (f> ^ Tp <;=4> (/, t) \= (f> if and only if {I , t) \= ip. 

The interpretation {I, r) is a model ofT, i.e., {I, r) |= T, iff {I, r) |= (j)for all (j) €T. 

'* We use typically symbols x, y, z to denote such free (functional) constants and symbols o, 6, c to denote proposi- 
tional atoms. 

^ The constants and operators appearing in a well-formed term t are based on a fixed width m. Moreover, the width 
w(a;) of each constant x € FC(T) must be the same throughout T. 



It is clear by Definition 2 that pure prepositional theories T are treated classically, i.e., (/, r) j= T iff 
/ 1= T in the sense of prepositional logic. As regards the theory T from Example 4, we have the sets of 
symbols At(T) = {a, b} and FC(T) = {x, y}. Furthermore, we observe that there is no model of T of 
the form {{a, b}, r) because it is impossible to satisfy x <2y and y <2 x simultaneously using any partial 
function r. On the other hand, there are 6 models of the form {{a}, r) because x <2 y can be satisfied in 
3 + 2 + 1 = 6 ways by picking different values for the 2-bit vectors x and y. 

3 Translation 

In this section, we present a translation of a logic program P into a bit- vector theory BV(P) that is similar 
to an existing translation [20] into difference logic. As its predecessor, the translation BV(P) consists of 
two parts. Clark's completion [8], denoted by CC(P), forms the first part of BV(P). The second part, i.e., 
R(P), is based on ranking constraints from [25] so that BV(P) = CC(P) U R(P). Intuitively, the idea is 
that the completion CC(P) captures supported models of P [1] and the further formulas in R(P) exclude 
the non-stable ones so that any classical model of BV(P) corresponds to a stable model of P. 
The completion CC(P) is formed for each atom a E At(P) on the basis of (2): 

1. If Defp(a) = 0, the formula is included to capture the corresponding empty disjunction in (2). 

2. If there is r e Defp(a) such that bd(r) = 0, then one of the disjuncts in (2) is trivially true and the 
formula a can be used as such to capture the definition of a. 

3. If Defp(a) = {r} for a rule r ^ P with n + m > 0, then we simplify (2) to a formula of the form 

/\ b A /\ ^c. (3) 

6ebd+(r) cebd-(r) 

4. Otherwise, the set Defp(a) contains at least two rules (1) with n + m > and 

y hdr (4) 

reDefp(a) 

is introduced using a new atom bd^. for each r G Defp(a) together with a formula 

bdr -H- /\ b A /\ ^c. (5) 

6ebd+(r) cebd-(r) 

The rest of the translation exploits the SCCs of the positive dependency graph of P that was defined in 
Section 2.1. The motivation is to limit the scope of ranking constraints which favors the length of the 
resulting translation. In particular, singleton components SCC(a) = {a} require no special treatment if 
tautological rules with a. e {bi. .... 6„} in (1) have been removed. Plain completion (2) is sufficient for 
atoms involved in such components. However, for each atom a e At(P) having a non-trivial component 
SCC(a) in DG'''(P) such that |SCC(a)| > 1, two new atoms exta and into are introduced to formahze 
the external and internal support for a, respectively. These atoms are defined in terms of equivalences 

V bd, (6) 

reExt_p(a) 

inta ^ V [b^^ ^ A ^°)] 

relntp(a) 66bd+(r)nSCC(a) 

where Xa and Xb are bit vectors of width m = [log2(|SCC(a)| + 1)] introduced for all atoms involved in 
SCC(a). The formulas (6) and (7) are called weak ranking constraints and they are accompanied by 

a exta V inta, (8) 
-.exta V -linta- (9) 

Moreover, when Ext p (a) and the atom a happens to gain external support from these rules, the value 
of Xa is fixed to by including the formula 



exta {Xa =m 0). 



(10) 



Example 5. Recall the program P from Example 1. The completion CC(P) is: 

a bdi V bd2. bdi^6Ac. bd2 o d. 

b -H- bda V bd4. bda ^ aA^d. bd4 o a A -ic. 

c o -id. 

Since P has only one non-trivial SCC, i.e., the component SCC(a) = SCC(6) = {a, b}, the weak ranking 
constraints resulting in R(P) are 

exto <-> bd2. inta -H- bdi A {xb <2 Xa). 
extfc -H- _L . 

intfc ^ [bds A (a;a <2 Xb)] V [bd4 A (xa <2 aj^)]. 

In addition to these, the formulas 

a — >■ exta V into, -lexta V -linta- exto {xa =2 0). 
b — > extb V iiitfa. -lextfc V -lint;,. 

are also included in R(_P). ■ 

Weak ranking constraints are sufficient whenever the goal is to compute only one answer set, or to check 
the existence of answer sets. However, they do not guarantee a one-to-one correspondence between the 
elements of AS(P) and the set of models obtained for the translation BV(P). To address this discrepancy, 
and to potentially make the computation of all answer sets or counting the number of answer sets more 
effective, strong ranking constraints can be imported from [20] as well. Actually, there are two mutually 
compatible variants of strong ranking constraints: 

bdr-)- V ^{Xb +mT^ <m Xa) (H) 

6ebd+(r)nSCC(a) 

inta^ V [bdr A V (Xa =m Xb +mT^)]- (12) 

relntp(a) bebd+(r)nSCC(a) 

The local strong ranking constraint (1 1) is introduced for each r € Intp(a). It is worth pointing out that 
the condition +m 1 <m Xa) is equivalent to Xb +m 1 >m Xa-^ On the other hand, the global variant 
(12) covers the internal support of a entirely. Finally, in order to prune copies of models of the translation 
that would correspond to the exactly same answer set of the original program, a formula 

^ {Xa =m 0) (13) 

is included for every atom a involved in a non-trivial SCC. We write R'(-P) and RS(P) for the respective 
extensions of R(P) with local/global strong ranking constraints, and R's(P) obtained using both. Similar 
conventions are applied to BV(P) to distinguish four variants in total. The correctness of these translations 
is addressed next. 

Theorem 1. Let P be a normal program and BV(P) its bit-vector translation. 

1. If S is an answer set of P, then there is a model {M, r) o/BV(P) such that S = M n At(P). 

2. If {M, t) is a model ofBV{P), then S = M n At(P) is an answer set of P. 

Proof. To estabUsh the correspondence of answer sets and models as formaUzed above, we appeal to the 
analogous property of the translation of P into difference logic (DL), denoted here by DL(P). In DL, theory 
atoms X < y + k constrain the difference of two integer variables x and y. Models can be represented as 
pairs (/, r) where 7 is a propositional interpretation and r maps constants of theory atoms to integers so 
that (/, r) \= X < y + k <^=^ t{x) < T{y) -\- k. The rest is analogous to Definition 2. 

* However, the form in (11) is used in our implementation, since -|-m and <m are amongst the base operators of the 
BOOLECTOR System. 



( ==> ) Suppose that S is an answer set of P. Then the results of [20] imply that there is a model 
(M, r) of DL(P) such that S = M D At{P). The valuation r is condensed for each non-trivial SCC S of 
DG~'~(P) as follows. Let us partition S into SqU . . . USn such that (i) T{xa) = t{xi,) for each < z < n 
and a,b £ Si, (ii) T{xa) = t{z)'^ for each a G Sq, and (iii) for each Q < i < j < n, a & Si, and h € Sj, 
T{xa) < T(a;6). Then define r' for the bit vector a;a associated with an atom a e by setting r'(a;a,i) = 1 
iff the j*^ bit of i is 1, i.e., T'{xa) = i- It follows that (J, r) \= Xb < Xa — 1 iff (/, t') \= Xb <m Xa for 
any a,b & S. Moreover, we have (M, r) ^ {xa < 2; + 0) A (2; < a;a + 0) iff (M, t') j= Xa =m for any 
a G 5. Due to the similar structures of DL(P) and BV(P), we obtain {M, r) |= BV(P) as desired. 

( ^ ) Let {M, t) be a model of BV(P). Then define r' such that t'{x) = Y.7=i{t{x, i) ■ 2^(^)-*) 
where x on the left hand side stands for the integer variable corresponding to the bit vector x on the right 
hand side. It follows that (/, r) |= Xb <m Xa iff r') \= Xb < Xa — I. By setting t'{z) = 0, we obtain 
(Af, r) 1= Xa =m if and only if {M, t') \= {xa < z + 0) A (z < Xa + 0). The strong analogy present in 
the structures of BV(P) and DL(P) imphes that (M, t') is a model of DL(P). Thus, S = Mn At{P) is 
an answer set of P by [20]. □ 

Even tighter relationships of answer sets and models can be established for the translations BV'(P), 
BVs(P), and BV's(P). It can be shown that the model {M, r) of BV*(P) corresponding to an answer set 
S' of P is unique, i.e., there is no other model {N, t') of the translation such that S = N Ci At(P). These 
results contrast with [20]: the analogous extensions DL*(P) guarantee the uniqueness of M in a model 
(M, r) but there are always infinitely many copies (M, r') of (M, r) such that (M, r') |= DL*(P). Such 
a valuation r' can be simply obtained by setting t'{x) = t{x) + 1 for any x. 



4 Native Support for Extended Rule Types 

The input syntax of the smodels system was soon extended by further rule types [28]. In solver interfaces, 
the rule types usually take the following simple syntactic forms: 

{ai, . . . ,ai} -s- 61, . . . ,bn, -ci, . . . ,~c„. (14) 
a Z{&i, . . . ,&„,~ci, . . . ,~c,„}. (15) 
a-(-l{bi = Wbi,...,b„ = Wb^ , ~ci = , . . . ,~Cm = Wc^}. (16) 

The body of a choice rule (14) is interpreted in the same way as that of a normal rule (1). The head, in 
contrast, allows to derive any subset of atoms ai, . . . ,ai, if the body is satisfied, and to make a choice 
in this way. The head a of a cardinality rule (15) is derived, if its body is satisfied, i.e., the number of 
satisfied literals amongst &i, . . . ,&„ and ^ci^ . . . ,^c„i is at least I acting as the lower bound. A weight 
rule of the form (16) generalizes this idea by assigning arbitrary positive weights to literals (rather than 
Is). The body is satisfied if the sum of weights assigned to satisfied literals is at least thus enabling one 
to infer the head a using the rule. In practise, the grounding components used in ASP systems allow for 
more versatile use of cardinality and weight rules, but the primitive forms (14), (15), and (16) provide a 
solid basis for efficient implementation via translations. The reader is referred to [28] for a generalization 
of answer sets for programs involving such extended rule types. The respective class of weight constraint 
programs (WCPs) is typically supported by SMODELS compatible systems. 

Whenever appropriate, it is possible to translate extended rule types as introduced above back to normjil 
rules. To this end, a number of transformations are addressed in [19] and they have been implemented as a 
tool called LP2normal'^. For instance, the head of a choice rule (14) can be captured in terms of rules 

ai <r- b, ~ai. ... ai <^ b, ^al. 
oT ^ai. ... a[-«— ~a;. 

where oT, ... ,07 are new atoms and 6 is a new atom standing for the body of (14) which can be defined 
using (14) with the head replaced by b. We assume that this transformation is applied at first to remove 

^ A special variable z is used as a placeholder for the constant in the translation DL(P) [20]. 
* http : // www . tcs . hut . fi/ Software/ asptools/ 



gringo program. Ip instance. Ip \ 

I smodels -internal -nolookahead \ 
I Ipcat -s=symbols . txt \ 
I lp2bv [-1] [-g] \ 
I boolector -fm 



Fig. 1. Unix shell pipeline for running a benchmark instance 



choice rules when the goal is to translate extended rule types into bit-vector logic. The strength of this 
transformation is locality, i.e., it can be applied on a rule-by-rule basis, and linearity with respect to the 
length of the original rule (14). To the contrary, linear normalization of cardinahty and weight rules seems 
impossible. Thus, we also provide direct translations into formulas of bit- vector logic. 

We present the translation of a weight rule (16) whereas the translation of a cardinality rule (15) is 
obtained as a special case Wh^ = . . . =wi,^ = Wci= ■ ■ ■ =Wc^ = 1. The body of a weight rule can be 
evaluated using bit vectors si, . . . of width k = \^og2^yll^=i + J2^i ''^a + 1)1 constrained by 

2 X {n + m) formulas 

bi (si =k Whl), -161 (si =fe 0), 

62 ->■ (S2 =fe Si +k Wb^), (S2 =k Si), 

K {Sn =k Sn-1 +k W^), -^K {Sn =k Sn-l), 

Ci ->■ {Sn+1 =k Sn), -'Cl (s„+i =fe S„ +k W^), 

^ {Sn+m ~k Sn+m—\)i ~^^'m ^ {.Sn+m ~k Sn+m—l ~t~fe ^Crn)* 

The lower bound I of (16) can be checked in terms of the formula -^{sn+m <fc where we assume that I 
is of width k, since the rule can be safely deleted otherwise. In view of the overall translation, the formula 
bdr <-> -i(s„+m <k 1) can be used in conjunction with the completion formula (4). Weight rules also 
contribute to the dependency graph DG~'~(P) in analogy to normal rules, i.e., the head a depends on all 
positive body atoms 61 , . . . ,6„. In this way, BV(P) generaUzes for programs P having extended rules. 



5 Experimental Results 

A new translator called lp2bv was implemented as a derivative of lp2diff^ that translates logic programs 
into difference logic. In contrast, the new translator will provide its output in the bit-vector format. In 
analogy to its predecessor, it expects to receive its input in the smodels file format. Models of the 
resulting bit-vector theory are searched for using boolector^' (v. 1.4.1) [6] and z3^^ (v. 2.11) [9] as 
back-end solvers. The goal of our preliminary experiments was to see how the performances of systems 
based on lp2bv compare with the performance of a state-of-the-art ASP solver clasp'^ (v. 1.3.5) [13]. The 
experiments were based on the NP-complete benchmarks of the ASP Competition 2009. In this benchmark 
collection, there are 23 benchmark problems with 516 instances in total. Before invoking a translator and 
the respective SMT solver, we performed a few preprocessing steps, as detailed in Figure 1, by calling: 

- GRINGO (v. 2.0.5), for grounding the problem encoding and a given instance; 

- smodels'^ (v. 2.34), for simplifying the resulting ground program; 

' http : // www . tcs . hut . f i/ Sof tware/lp2dif f / 
http : / /www . tcs . hut . f i/Sof tware/smodels/ 
" http : / /fmv . jku. at /boolector/ 

'^http://research. microsoft. com/ en-us/um/ redmond/pro ject s/ z3/ 
http : / / www . cs . uni-potsdam. de/ clasp/ 
http : / /www . tcs . hut . fi/ Software/ smodels/ 



Table 1. Experimental results without normalization 
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CLASP 




LP2BV+B00LECT0R 




LP2BV+Z3 






LP2DIFF+Z3 




Benchmark 










W 


L 


G 


LG 


w 


L 


G 


LG 


w 


L 


G 


LG 




































Overall Performance 


516 




465 




276 


244 


261 


256 


217 


216 


194 


204 


360 


349 


324 


324 








347/118 




188/88 


161/83 


174/87 


176/ 80 


142/75 


147/ 69 


124/70 


135/ 69 


257/103 


251/98 


225/ 99 


226/ 98 




































KnightTour 


10 




8/0 




2/0 


1/0 


0/0 


0/0 


1/0 


0/0 


0/0 


1/0 


6/0 


6/0 


4/0 


5/0 


GraphColouring 


29 




8/0 




7/0 


7/0 


7/0 


7/0 


6/0 


7/0 


7/0 


7/0 


7/0 


7/0 


7/0 


7/0 


WireRouting 


23 




11/11 




2/3 


1/ 1 


1/2 


0/2 


1/3 


0/0 


0/0 


0/ 1 


3/3 


2/3 


2/4 


5/3 


DisjunctiveScheduling 


10 




5/0 




0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


GraphPartitioning 


13 




611 




3/0 


3/0 


3/0 


3/0 


4/0 


4/0 


4/0 


3/0 


6/2 


6/ 1 


6/ 1 


6/ 1 


ChannelRouting 


11 




6/2 




6/2 


6/2 


6/2 


6/2 


5/2 


6/2 


6/2 


6/2 


6/2 


6/2 


6/2 


6/2 


Solitaire 


27 




19/0 




2/0 


5/0 


1/0 


4/0 


0/0 


0/0 


0/0 


0/0 


21/0 


21/0 


20/0 


21/0 


Labyrinth 


29 




26/0 




1/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


WeightBoundedDominatingSet 


29 




26/0 




18/0 


18/0 


17/0 


18/0 


12/0 


12/0 


11/0 


12/0 


22/0 


22/0 


22/0 


21/0 


MazeGeneration 


29 




10/15 




8/15 


1/15 


0/15 


0/16 


5/16 


1/15 


0/15 


1/15 


10/17 


10/15 


5/15 


4/15 


ISPuzzle 


16 




16/0 




16/0 


15/0 


14/0 


15/0 


4/0 


4/0 


5/0 


5/0 


0/0 


0/0 


0/0 


0/0 


BlockedNQueens 


29 




15/14 




2/2 


0/2 


1/2 


0/2 


1/0 


2/0 


2/0 


0/0 


15/13 


15/13 


15/12 


15/13 


ConnectedDominatingSet 


21 




10/10 




10/11 


9/8 


10/11 


6/3 


10/10 


9/10 


10/9 


10/9 


9/8 


7/6 


9/7 


7/6 


EdgeMatching 


29 




29/0 




0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


0/0 


3/0 


1/0 


3/0 


2/0 


Fastfood 


29 




10/19 




9/16 


10/16 


10/16 


9/16 


9/9 


9/9 


9/10 


9/9 


10/18 


10/18 


10/18 


10/18 


GeneralizedSlitherlink 


29 




29/0 




29/0 


20/0 


29/0 


29/0 


29/0 


29/0 


16/0 


29/0 


29/0 


29/0 


29/0 


29/0 


HamiltonianPath 


29 




29/0 




27/0 


25/0 


29/0 


28/0 


26/0 


27/0 


25/0 


26/0 


29/0 


29/0 


29/0 


29/0 


Hanoi 


15 




15/0 




15/0 


15/0 


15/0 


15/0 


5/0 


5/0 


5/0 


4/0 


15/0 


15/0 


15/0 


15/0 


HierarchicalClustering 


12 




8/4 




8/4 


8/4 


8/4 


8/4 


4/4 


4/4 


4/4 


4/4 


8/4 


8/4 


8/4 


8/4 


SchurNumbers 


29 




13/16 




6/16 


5/16 


5/16 


5/16 


9/16 


9/16 


9/16 


9/16 


11/16 


11/16 


11/16 


11/16 


Sokoban 


29 




9/20 




9/19 


8/19 


8/19 


8/19 


7/15 


7/13 


7/14 


5/13 


9/20 


9/20 


9/20 


9/20 


Sudoku 


10 




10/0 




5/0 


4/0 


4/0 


5/0 


4/0 


4/0 


4/0 


4/0 


9/0 


8/0 


8/0 


9/0 


TravellingSalesperson 


29 




29/0 




3/0 


0/0 


6/0 


10/0 


0/0 


8/0 


0/0 


0/0 


29/0 


29/0 


7/0 


7/0 



- LPCAT (v. 1.18), for removing all unused atom numbers, for making the atom table of the ground 
program contiguous, and for extracting the symbols for later use; and 

- LP2normal (version 1.11), for normalizing the program. 

The last step is optional and not included as part of the pipeline in Figure 1. Pipelines of this kind were 
executed under Linux/Ubuntu operating system miming on six-core AMD Opteron^'^'^) 2435 processors 
under 2.6 GHz clock rate and with 2.7 GB memory linait that corresponds to the amount of memory 
available in the ASP Competition 2009. 

For each system based on a translator and a back-end solver, there are four variants of the system to con- 
sider: W indicates that only weak ranking constraints are used, while L, G, and LG mean that either local, 
or global, or both local and global strong ranking constraints, respectively, are employed when translating 
the logic program. 

Table 1 collects the results from our experiments without normalization whereas Table 2 shows the 
results when lp2normal [19] was used to remove extended rule types discussed in Section 4. In both 
tables, the first colunm gives the name of the benchmark, followed by the number of instances of that 
particular benchmark in the second column. The following columns indicate the numbers of instances that 
were solved by the systems considered in our experiments. A notation like 8/4 means that the system was 
able to solve eight satisfiable and four unsatisfiable instances in that particular benchmark. Hence, if there 
are 15 instances in a benchmark and the system could only solve 8/4, this means that the system was 
unable to solve the remaining three instances within the time limit of 600 seconds, i.e. ten minutes, per 
instance'^. As regards the number of solved instances in each benchmark, the best performing translation- 
based approaches are highlighted in boldface. Though it was not shown in all tables, we also run the 
experiments using translator LP2DIFF with z3 as back-end solver, and the summary is included in Table 
3 — giving an overview of experimental results in terms of total numbers of instances solved out of 516. 

It is apparent that the systems based on LP2bv did not perform very well without normalization. As 
indicated by Table 3, the overall performance was even worse than that of systems using LP2DIFF for 

One observation is that the performance of systems based on lp2bv is quite stable: even when we extended the 
time limit to 20 minutes, the results did not change much (differences of only one or two instances were perceived 
in most cases). 



Table 2. Experimental results with normalization 





INST 


CLASP 


LP2BV+BOOLECTOR 




LP2BV+Z3 




Benchmark 






W 


L 


G 


LG 


W 


L 


G 


LG 


Overall Performance 


516 


459 


381 


343 


379 


381 


346 


330 


325 


331 






346/113 


279/102 


243/100 


278/101 


281/100 


240/106 


231/99 


224/101 


232/ 99 


KnightTour 


10 


10/0 


2/0 


2/0 


1/0 


0/0 


1/0 


0/0 


0/0 


0/0 


GraphColouring 


29 


9/0 


8/0 


8/ 


8/0 


8/0 


9/2 


9/2 


9/2 


9/2 


WireRouting 


23 


11/11 


2/6 


1/3 


1/3 


1/3 


2/7 


1/4 


1/4 


1/3 


DisjunctiveScheduling 


10 


5/0 


S/0 


5/0 


5/0 


5/0 


5/0 


5/0 


5/0 


5/0 


GraphPartitioning 


13 


4/ 1 


5/0 


5/0 


4/0 


5/0 


2/ 1 


2/ 1 


2/ 1 


2/0 


ChannelRouting 


11 


6/2 


6/2 


6/2 


6/2 


6/2 


6/2 


6/2 


6/2 


6/2 


Solitaire 


27 


18/0 


23/0 


23/0 


23/0 


23/0 


22/0 


22/0 


22/0 


22/0 


Labyrinth 


29 


27/0 


1/0 


1/0 


2/0 


3/0 


0/0 


0/0 


0/0 


0/0 


WeightBoundedDominatingSet 


29 


25/0 


15/0 


15/0 


15/0 


16/0 


10/0 


10/0 


10/0 


10/0 


MazeGeneration 


29 


10/15 


8/15 


0/15 


0/15 


0/16 


5/16 


0/15 


0/15 


0/15 


15Puzzle 


16 


15/0 


16/0 


16/0 


16/0 


16/0 


11/0 


10/0 


11/0 


11/0 


BlockedNQueens 


29 


15/14 


14/14 


14/14 


14/14 


14/14 


15/14 


15/14 


15/14 


15/14 


ConnectedDominatingSet 


21 


10/11 


10/11 


8/11 


9/11 


9/10 


10/11 


9/11 


9/11 


9/11 


EdgeMatching 


29 


29/0 


29/0 


29/0 


29/0 


29/0 


29/0 


29/0 


29/0 


29/0 


Fastfood 


29 


10/19 


9/14 


9/15 


9/16 


9/15 


0/13 


0/10 


0/12 


0/12 


GeneralizedSlitherlink 


29 


29/0 


29/0 


21/0 


29/0 


29/0 


29/0 


29/0 


21/0 


29/0 


HamiltonianPath 


29 


29/0 


29/0 


28/0 


29/0 


29/0 


29/0 


29/0 


29/0 


29/0 


Hanoi 


15 


15/ 


15/0 


15/0 


15/0 


15/0 


15/0 


15/0 


15/0 


15/0 


HierarchicalClustering 


12 


8/4 


8/4 


8/4 


8/4 


8/4 


8/4 


8/4 


8/4 


8/4 


SchurNumbers 


29 


13/16 


10/16 


10/16 


9/16 


10/16 


13/16 


13/16 


13/16 


13/16 


Sokoban 


29 


9/20 


9/20 


9/20 


9/20 


9/20 


9/20 


9/20 


9/20 


9/20 


Sudoku 


10 


10/0 


10/0 


10/0 


10/0 


10/0 


10/0 


10/0 


10/0 


10/0 


TravellingSalesperson 


29 


29/0 


16/0 


0/0 


27/0 


27/0 


0/0 


0/0 


0/0 


0/0 



translation and z3 for model search. However, if the input was first translated into a normal logic program 
using lp2normal, i.e., before translation into a bit-vector theory, the performance was clearly better. 
Actually, it exceeded that of the systems based on lp2diff and became closer to that of CLASP. We note 
that normalization does not help so much in case of lp2diff and the experimental results obtained using 
both normalized and unnormalized instances are quite similar in terms of solved instances. Thus it seems 
that solvers for bit- vector logic are not able to make the best of native translations of cardinality and weight 
rules from Section 4 in full. If an analogous translation into difference logic is used, as implemented 
in lp2diff, such a negative effect was not perceived using z3. Our understanding is that the efficient 
graph-theoretic satisfiability check for difference constraints used in the search procedure of z3 turns the 
native translation feasible as well. As indicated by our test results, BOOLECTOR is clearly better back-end 
solver for lp2bv than z3. This was to be expected since boolector is a native solver for bit-vector 
logic whereas z3 supports a wider variety of SMT fragments and can be used for more general purposes. 
Moreover, the design of LP2bv takes into account operators of bit- vector logic which are directly supported 
by BOOLECTOR and not implemented as syntactic sugar. 

In addition, we note on the basis of our results that the performance of the state-of-the-art ASP solver 
CLASP is significantly better, and the translation-based approaches to computing stable models are still 
left behind. By the results of Table 2, even the best variants of systems based on LP2bv did not work well 
enough to compete with clasp. The difference is especially due to the following benchmarks: Knight Tour, 
Wire Routing, Graph Partitioning, Labyrinth, Weight Bounded Dominating Set, Fastfood, and Travelling 
Salesperson. All of them involve either recursive rules (Knight Tour, Wire Routing, and Labyrinth), weight 
rules (Weight Bounded Dominating Set and Fastfood), or both (Graph Partitioning and Travelling Sales- 
person). Hence, it seems that handling recursive rules and weight constraints in the translational approach 
is less efficient compared to their native implementation in clasp. When using the current normalization 
techniques to remove cardinality and weight rules, the sizes of ground programs tend to increase signifi- 
cantly and, in particular, if weight rules are abundant. For example, after normalization the ground programs 
are ten times larger for the benchmark Weight Bounded Dominating Set, and five times larger for Fastfood. 
It is also worth pointing out that the efficiency of CLASP turned out to be insensitive to normaUzation. 



Table 3. Summary of the experimental results 



System 


W 


L 


G 


LG 


LP2BV+BOOLECTOR 


276 


244 


261 


256 


LP2BV+Z3 


217 


216 


194 


204 


LP2DIFF+Z3 


360 


349 


324 


324 


CLASP 


465 


LP2NORMAL2BV+BOOLECTOR 


381 


343 


379 


381 


LP2NORMAL2BV+Z3 


346 


330 


325 


331 


LP2NORMAL2D1M+Z3 


364 


357 


349 


349 


LP2NORMAL+CLASP 


459 



While having trouble with recursive rules and weight constraints for particular benchmarks, the transla- 
tional approach handles certain large instances quite well. The largest instances in the experiments belong 

to the Disjunctive Scheduling benchmark, of which all instances are ground programs of size over one 
megabyte but after normalization'^, the lp2bv systems can solve as many instances as clasp. 

6 Conclusion 

In this paper, we present a novel and concise translation from normal logic programs into fixed-width bit- 
vector theories. Moreover, the extended rule types supported by smodels compatible answer set solvers 
can be covered via native translations. The length of the resulting translation is Unear with respect to the 
length of the original program. The translation has been implemented as a translator, LP2bv, which enables 
the use of bit-vector solvers in the search for answer sets. Our preliminary experimental results indicate 
a level of performance which is similar to that obtained using solvers for difference logic. However, this 
presumes one first to translate extended rule types into normal rules and then to apply the translation into 
bit- vector logic. One potential explanation for such behavior is the way in which SMT solvers implement 
reasoning with bit vectors: a predominant strategy is to translate theory atoms involving bit vectors into 
propositional formulas and to apply satisfiability checking techniques systematically. We anticipate that an 
improved performance could be obtained if a native support for certain bit vector primitives were incor- 
porated into SMT solvers directly. When comparing to the state-of-the-art ASP solver clasp, we noticed 
that the performance of the translation based approach compared unfavorably, in particular, for benchmarks 
which contained recursive rules or weight constraints or both. This indicates that the performance can be 
improved by developing new translation techniques for these two features. In order to obtain a more com- 
prehensive view of the performance characteristics of the translational approach, the plan is to extend our 
experimental setup to include benchmarks that were used in the third ASP competition [7]. Moreover, we 
intend to use the new SMT library format [4] in future versions of our translators. 
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